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Abstract. Despite extensive work on the interplay between traffic dynamics and 
epidemic spreading, the control of epidemic spreading by routing strategies has not 
received adequate attention. In this paper, we study the impact of efficient routing 
protocol on epidemic spreading. In the case of infinite node-delivery capacity, where 
the traffic is free of congestion, we find that that there exists optimal values of routing 
parameter, leading to the maximal epidemic threshold. This means that epidemic 
spreading can be effectively controlled by fine tuning the routing scheme. Moreover, 
we find that an increase in the average network connectivity and the emergence of 
traffic congestion can suppress the epidemic outbreak. 
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1. Introduction 

Epidemic spreading [H|2J|3JIU|5J|6J[7J|HJ|9J ® HD H2] and traffic dynamics [131 HU H5l 
m im nm urn m] on complex networks have attracted much attention in the past 
decade. For a long time, the two types of dynamical processes have been studied 
independently. However, in many cases, epidemic spreading is relied on the process 
of transportation. For example, a computer virus can spread over Internet via data 
transmission [HI [22]. Another example is that air transport tremendously accelerates 
the propagation of infectious diseases among different countries. 

The first attempt to incorporate traffic into epidemic spreading is based on 
metapopulation model [23] EH EE, EH E3 EH EH [30] EH [32] . This framework describes 
a set of spatially structured interacting subpopulations as a network, whose links denote 
the traveling path of individuals across different subpopulations. Each subpopulation 
consists of a large number of individuals. An infected individual can infect other 
individuals in the same subpopulation. In a recent work, Meloni et al. proposed another 
traffic-driven epidemic spreading model [33], in which each node of a network represents 
a router in the Internet and the epidemic can spread between nodes by the transmission 
of packets. A susceptible node will be infected with some probability every time it 
receives a packet from an infected neighboring node. 

Despite broad interests in traffic-driven epidemic spreading, the control of epidemic 
spreading by routing strategies has received little attention. In another recent work, 
Meloni et al. observed that when travelers decide to avoid locations with high levels 
of prevalence, this self-initiated behavioral change may enhance disease spreading [33]. 
Later, Yang et al. found that epidemic spreading can be effectively controlled by a local 
routing strategy [35]. In the local routing protocol, each node does not know the whole 
network’s topological information and the packet is forwarded to a neighboring node i 
with a probability that is proportional to the power of Ts degree [36]. It is noted that 
in the local traffic routing, the average traveling time of a packet (T) is proportional to 
the network size N [37]. However, in global routing protocols such as the shortest-path 
routing, (T) usually increases approximately logarithmically with N [38]. Thus, from 
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the view of transmission time, global routing protocols may be superior to local routing 
protocols. 

So far, the control of traffic-driven epidemic spreading by a global routing protocol 
has not been studied. To address the above issue, we consider an efficient routing 
strategy proposed by Yan et al. [39]. In the efficient routing protocol, each node in a 
network is assigned a weight that is proportional to the power of its degree, where the 
power exponent a is a tunable parameter. The efficient path between any two nodes is 
corresponding to the route that makes the sum of the nodes’ weight (along the path) 
minimal. It has been proved that the traffic throughput of the network can be greatly 
improved by employing the efficient routing strategy as compared to the shortest-path 
strategy [39]. In this paper, we intend to study how the the efficient routing protocol 
affects traffic-driven epidemic spreading. Our preliminary results have shown that there 
exists an optimal value of a, leading to the maximal epidemic threshold. 

The paper is organized as follows. In Sec. [2j we formalize the problem by 
introducing the efficient routing strategy into traffic-driven epidemic spreading. In 
Sec. [3] we investigate the epidemic spreading on scale-free networks by considering 
two cases of node-delivering capacity, i.e., infinite capacity and finite capacity. Finally, 
conclusions and discussions are presented in Sec. [H 

2. Model and Methods 

Following the work of Meloni et al. [33], we incorporate the traffic dynamics into the 
classical susceptible-infected-susceptible model of epidemic spreading as follows. 

(i) Efficient routing protocol. In a network of size N, at each time step, A N new 
packets are generated with randomly chosen sources and destinations (we call A as the 
packet-generation rate), and each node can deliver at most C packets towards their 
destinations. For any path between nodes i and j, P{i — > j) := i = X\,X 2 ,x n = j, 
we define 

n 

L{P{i->j):a) = ^2k(xi) a , ( 1 ) 

i=i 

where k{x{) is the degree of node xi and a is a tunable parameter. For any given a, 
the efficient path between i and j is corresponding to the route that makes the sum 
L ( P(i — > j ) : a) minimum. Packets are delivered following the efficient path. When 
a = 0, the efficient path recovers the traditional shortest path. Once a packet reaches 
its destination, it is removed from the system. The queue length of each node is assumed 
to be unlimited and the first-in-first-out principle holds for the queue. 

(ii) Epidemic dynamics. After a transient time, the total number of delivered 
packets at each time will reach a steady value, then an initial fraction of nodes p 0 is set 
to be infected (e.g., we set po = 0.1 in numerical experiments). The infection spreads in 
the network through packet exchanges. Each susceptible node has the probability fd of 
being infected every time it receives a packet from an infected neighbor. The infected 
nodes recover at rate /i (we set pi — 1 in this paper). 
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Figure 1. The epidemic threshold /3c as a function of a for different values of ( k ). 
The symbols are from numerical simulations. The solid, dashed, and dotted curves 
correspond to the theoretical prediction from Eq. (2) for ( k ) = 4, 6 and 10, respectively. 
The inset shows the optimal value a op t of the routing parameter as a function of node 
degree ( k). The packet-generation rate A = 1. Each data point results from an average 
over 100 different realizations. 


3. Main results and Analysis 

In the following, we carry out simulations systematically by employing traffic-driven 
epidemic spreading on the Barabasi-Albert (BA) scale-free networks [4l]J. The size of 
the BA network is set to be N = 2000. In the case where the node-delivering capacity 
is infinite (C —> oo), traffic congestion will not occur in the network. When the node- 
delivering capacity is finite, traffic congestion can occur if the packet-generating rate A 
exceeds a critical value [35]. Therefore infinite and finite node-delivering capacity are 
considered respectively in the following sections. 

3.1. Infinite node-delivering capacity 

Previous studies have shown that there exists an epidemic threshold (3 C , below which the 
epidemic goes extinct [33, 135]. Figure Q] shows the dependence of /3 C on a for different 
values of the average degree ( k ) of the network. We fold that for each value of (k), 
there exists an optimal value of a, hereafter denoted by a opt , leading to the maximum 
/ 3 C . The inset of Fig. |T] shows that a op t decreases from 0.7 to 0.5 as ( k ) increases from 4 
to 16. 

According to the analysis in Ref. [33], the epidemic threshold for uncorrelated 
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Figure 2. The epidemic threshold /3 C as a function of (k) for different values of a. 
The packet-generation rate A = 1. Each data point results from an average over 100 
different realizations. 


networks is 

a _ (fraig) 1 
(^ig)A N’ 

where 6 a i g is the efficient algorithmic betweenness of a node p®, 33] and (•) denotes 
the average of all nodes. The efficient algorithmic betweenness of a node represents 
the average number of packets passing through that node at each time step when the 
packet-generation rate A = 1/iV. In this paper, the efficient algorithmic be- tweenness 
of a node k can be calculated as 

.fc _ 1 'V' C 

a tJ 

where is the total number of efficient paths going from i to j, and <Jij(k) is the 
number of efficient paths going from i to j and passing through k. From Fig. [H one can 
see that the theoretical predictions agree with numerical results qualitatively. 

Next, we study the effect of the average degree of the network on the traffic-driven 
epidemic spreading. Figure [2] shows the epidemic threshold j3 c as a function of the 
average degree of the network ( k) for different values of a. From Fig. [2J we fold that 
for each value of a, /3 C increases with ( k ), in contrast to the behavior of spreading 
dynamics in the absence of traffic pH]. This phenomenon can be understood as follows. 
An increase in the average degree of the network shortens the average time steps that a 
packet spends traveling from its source to its destination and decreases the number of 
packages passing through each node, leading to a decrease in the infection probability 
of each node. 

Figure [3] shows the epidemic threshold (3 C as a function of the packet-generation rate 
A for different values of a. One can see that /3 C scales inversely with A, as predicted by 
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Figure 3. The epidemic threshold j3 c as a function of the packet-generation rate A 
for different values of a. The average degree of the network (k) = 6. The slope of 
the fitted line is about -1. Each data point results from an average over 100 different 
realizations. 



Figure 4. The infection probability of nodes p as a function of the algorithmic 
betweenness 6 a i g for different values of a. The average degree of the network (k) = 6. 
Each data point results from an average over 100 different realizations. 


Eq. (2), indicating that the increase of traffic flow facilitates the outbreak of epidemic. 
The similar result has also been found in Ref. |33j . 

An interesting issue is how the algorithmic betweenness & a i g affects the infection 
probability of nodes p. Figure [4] features the dependence of p on 6 alg for different values 
of a. From Fig. [2 one can observe that the algorithmic betweenness of the nodes is 
positively correlated with the risk of them being infected. 
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Figure 5. The epidemic threshold /3 C as a function of the packet-generation rate A 
for finite and infinite C. The average degree of the network ( k) = 6 and the routing 
parameter a = 0.7. The critical packet-generating rate for C = 10 is A c « 0.7. Each 
data point results from an average over 100 different realizations. 


3.2. Finite node-delivering capacity 

In the case of finite node-delivering capacity, traffic congestion occurs when the packet¬ 
generating rate exceeds a critical value A c , which can can be estimated as [331135], 

C 

^ c ~ yyrmax ’ (^) 

alg 

where b^ x is the largest algorithmic betweenness of the network. 

Figure [5] shows the epidemic threshold /3 C as a function of the packet-generation rate 
A for finite and infinite C. One can see that when A < A c , j3 c is identical for both cases of 
the finite and infinite delivery capacities. However, for A > A c , /3 C is larger in the case of 
finite capacity than that in the infinite capacity case, indicating that traffic congestion 
can suppress the outbreak of epidemic. This result is consistent with that in Ref. [ 33] . 
The above phenomenon can be explained as follows. Once a node becomes congested, 
it cannot deliver the total packets in its queue at each time step. A decrease in the 
number of delivered packets can help nodes reduce the probability of being infected. 

Figure [6] shows the epidemic threshold /3 C as a function of a for different values of 
the packet-generation rate A. One can observe that for small values of A (e.g., A = 1 
or A = 1.5), there exists an optimal value of a, leading to the maximal (3 C . However, 
when A is large enough (e.g., A = 4), j3 c decreases with the increase of a. The inset of 
Fig. [6] shows the critical packet-generating rate A c as a function of a. One can see A c is 
maximized by an optimal value of a. 
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Figure 6. The epidemic threshold /3 C as a function of a for different values of the 
packet-generation rate A. The average degree of the network (k) = 6 and the node- 
delivering capacity C = 50. The inset shows that the critical packet-generating rate 
A c as a function of a. When A = 1, A > A c for a < 0.4 and a > 2.2. When A = 1.5, 
A > A c for a < 0.5 and a > 1.6. When A = 4, A > A c for all the values of a. Each 
data point results from an average over 100 different realizations. 


4. Conclusions and Discussions 

In conclusion, we have studied the impact of efficient routing protocol on traffic-driven 
epidemic spreading. We find that the epidemic threshold increases with the average 
degree of the network when other parameters are fixed. Besides, we find that nodes 
with larger algorithmic betweenness are more likely to be infected. Both analytic 
and numerical results show that, there exists an optimal value of routing parameter, 
leading to the maximal epidemic threshold. This means that epidemic spreading can 
be controlled by fine tuning the routing scheme. We hope our results can be useful to 
understand and control spreading dynamics. 
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